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, Abstract. Long Period Variable stars are ubiquitous among the bright red giant branch stars in NGC 5128. Mostly 

they are found to be brighter than the tip of the first ascent giant branch with magnitudes ranging from about 
K — 19 to K = 21.5. They have periods between 155 and 1000 days and K-band amplitudes between 0.1 and 2 
mag, characteristic of semi-regular and Mira variables. We compare the colors, periods and amplitudes of these 
^ ' variables with those found in old stellar populations like Galactic globular clusters and Galactic bulge as well as 

OO I with intermediate-age Magellanic Cloud long period variables. The population of stars above the tip of the red 

lO ■ giant branch (RGB) amounts to 2176 stars in the outer halo field (Field 1) and 6072 stars in the inner halo field 

■ (Field 2). The comparison of the luminosity functions of the Galactic bulge, M31 bulge and NGC 5128 halo fields 

shows an excess of bright AGB stars extending to Mk — —8.65. The large majority of these sources belong to the 
asymptotic giant branch (AGB) population in NGC 5128. Subtracting the foreground Galactic stars and probable 
blends, at least 26% and 70% of AGB stars are variable in Fields 1 and 2, respectively. The average period of 
NGC 5128 LPVs is 395 days and the average amplitude 0.77 mag. Many more short period Miras are present 
• in Field 2 than in Field 1 indicating a difference in the stellar populations between the two fields. Period and 

amplitude distributions and near-IR colors of the majority of LPVs in NGC 5128 are similar to the Galactic bulge 
Q . variables. However, some ~ 10% of LPVs have periods longer than 500 days and thus probably more massive, 

' hence younger, progenitor stars. A few carbon star candidates are identified based on their red J ~ H and H — K 

' colors. 

■ 
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1. Introduction & Ellas Frogel & Whitelock [TM5|l . The frequency of 

LPVs in old metal-rich globular clusters of the MW and in 

In intermediate-age populations (-1-5 Gyr old) numer- ^-^^ g^jg^ ^^^^ ^^^^-^^ p^^g^l ^ Whitelock fTTM . 

ous bright asymptotic giant branch (AGB) stars are lo- Qld populations of lower metallicity are known not to have 

cated above the tip of the RGB. However, bright stars have ^GB stars brighter than the RGB tip. 

also been detected above the tip of the RGB among old mi i • w _l r-mr • 1 1 ^ n i i-i 

, ,. „ , , , , , , •,, m /TTi ^ i he brightest LPVs m old stellar populations hke the 

populations hke Galactic globular clusters with be/H <, ^ , ,. , , , , ; , n^-n 

„ . , • ,, ^ , ,. , , 1 D T^i- -.r^oo Galactic globular clusters and the old disk m the Milky 

— 1.0 dex and m the Galactic bulge (Lrogel & Flias 1988, T ,^ o /t>«- • , i rmrrm . , 

• • X 1 rrmwi f i Vtwh, • i • Way reach Mk = -8 (Mennessier et al. 120011 Feast et 

Guarnieri et al. 119971 Momany et al. l2003ll . implying the , rwrm t.^ ^ ^ iaaaai , , • , /, 

ri • iTTTTti ^ • . , . , ' , , / , al. 120021 Momany et al. I2003II . The higher mass (hence 

presence ot bright AGB stars m metal-rich and old popula- \ x m r • . i t ■» i ■ , i i • i r i 

, p , . , , . , , ,, „^Ti . • • 11 younger) LPVs m the LMG and m the disk oi our galaxy 

tions. All or the bright giants above the RGB tip m globu- ° ' ^ , . , r ^ wtt , c 

, , , , 1 , -1 • , 1 1 can be more luminous, reaching Mk — —9.4 (Hughes & 

lar clusters seem to be long period variables (LPVs; Lrogel ,iaaa • ^ i riS,.i , 

° \ J o Wood 1990. Mennessier et al. I2D01) . 

'Z Z £f ~i , , A J r> -1 u There exists a long standing debate about the presence 

bend ojjprint requests to: M. Kejkuba o o r 

* Based on observations collected at the European of intermediate-age population in dwarf elliptical galaxies 

Southern Observatory, Paranal, Chile, within the Observing and spiral bulges. On the basis of the Wide Field Camera 

Programmes 63.N-0229, 65.N-0164, 67.B-0503, 68.B-0129 and 1 (WF Cl) H ST data, Holtzman et al. H1993|l and Vallenari 

69.B-0292, and at La Silla Observatory, Chile, within the et al. fl99Ci| concluded that the majority of stars in the 

Observing Programme 64.N-0176(B). Galactic bulge are of intermediate age. Other studies in- 



2 



Rejkuba et al.: LPVs in Cen A 



dicated mainly old metal-rich population (e.g. Frogel et 
al. 119901 Tiede et al. 11995(1 . Newer deep optical and near- 
IR CMDs suggest that the Bulge is old without even a 
trace of an intermediate-age population (e.g. Ortolani et 
al. [1995, Feltzing & Gilmore i^OOl Zoccah et al. pOS ). 
In the M32 dwarf elliptical galaxy the situation might be 
different, given the presence of an extended giant branch 
reaching Mk — —8.7 {Mboi — —5.5), indicative of an 
intermediate-age 4 Gyr) population, which has been 
revealed by near-IR photometry by Elston & Silva l|1992|l 
and Freedman H1992|l . However, it is of interest to note 
that in the deep optical VI photometry of M32 with HST, 
no optically bright AGB stars have been found (Grillmair 
et al. I1996|l . which might be expected due to large bolo- 
metric corrections. 

In the gE NGC 5128, Soria et al. ifTU^ . suggested 
the presence of up to 10% bright AGB stars belonging 
to an intermediate-age population, based on a VI HST 
CMD. Marleau et al. ((?000)1 made a similar suggestion 
based on JH NICMOS data. Harris et al. (1999, 2000), 
on the contrary, do not find any bright AGB stars in their 
VI CMDs of two halo fields in NGC 5128. As mentioned 
above, V and / bands are not very sensitive indicators of 
these cool giants and thus some of them might have been 
confused with first ascent red giants, foreground stars and 
few stellar blends or stars with larger photometric errors. 
In the two fields observed in V and ^fC-bands with VLT 
(Rejkuba et al. i20DIfl . a large number of bright AGB stars 
have been detected extending up to bolometric magni- 
tude of —5. Most of these stars have been found to be 
variables in our long term monitoring programme with 
ISAAC at VLT. The LPV catalogue has been published by 
Rejkuba et al. H2003|l . Here we analyse near-IR properties 
of these variables and compare them with the LPV popu- 
lation found in the Magellanic Clouds and in the Galactic 
bulge, with the aim of constraining the contribution of 
intermediate-age stars to the NGC 5128 halo population. 

2. The Catalogue 

The LPV variables catalogue in NGC 5128 has been com- 
piled and published by Rejkuba et al. (2003). The reader 
is referred to that paper for the details on the data re- 
duction and the determination of periods. In summary, 
the observations consist of single-epoch 1-hour exposures 
in Js and H-hands and multi-epoch photometry in Ks- 
band obtained with ISAAC near-IR array at ESQ Paranal 
Observatory with UTl VLT. The two fields observed are 
located 17' north-east (Field 1) and 9' south (Field 2) 
from the center of NGC 5128. 

i^Ts-band observations span a time interval of 1197 
days, from April 1999 till July 2002, with 20 epochs in 
Field 1 and 23 in Field 2. One additional 45 min observa- 
tion of Field 2 in if^-band with SOFI at the NTT at La 
Silla Observatory brings the total number of epochs to 24 
for that field. It should be noted that SOFI near-IR array 
is a scaled version of ISAAC and that the Kg-hand filters 
at the two instruments are identical. 



The photometry of this homogeneous data set has been 
performed with DAOPHOT and ALLFRAME software 
(Stetson Ll9Mj) . The 50% completeness limits in Kg and 
i?-bands are at 22.5 for Field 1 and 21.5 for Field 2. 
In Js-band 50% completeness is achieved at Js ~ 23.25 
for Field 1 and Js = 22.5 for Field 2. To aU X.-band 
magnitudes from the catalogue (Rejkuba et al. i2003j) we 
have subtracted 0.1 mag. This was necessary in order to 
match the observed colors of foreground dwarf and giant 
stars with that of Bessell & Brett (|1988|l fiducials (see 
also Sect. 17. H and Reikuba l2003|l . We have re-checked the 
photometric calibration and most probably this 0.1 offset 
in K-hand is due to an error in the aperture correction. 
It should be noted that the aperture corrections in such 
crowded fields are rather uncertain and very difficult to 
measure accurately. 

There are total of 15574 sources in Field 1 and 18098 
in Field 2 that were detected in at least 3 Kg-hand epochs. 
Restricting the detection to Js, H and 3 Ks-h&nd epochs 
the total number of objects is 13111 and 16434 in Fields 
1 and 2, respectively. Among these there are more than 
1500 variable stars. For 1046 red variables, with at least 
10 i^Tg-band measurements, periods and amplitudes were 
determined with Fourier analysis and then refined with 
non-linear fitting of sinusoidal function: 

K{t) = A cos (^2^ii^^ + B sin (^27r^^^^ + A'o(l) 

In the final catalogue (Rejkuba et al. I2003|l periods 
of variable stars are average values determined from the 
Fourier analysis and the non-linear sine-curve fitting algo- 
rithms, except in the cases where a visual inspection of the 
light curve clearly preferred one of the two cases. In most 
cases the two periods were nearly equal (see Rejkuba et 
al. 120031 for simulations and detailed discussion of period 
accuracy). Amplitudes, a = 2 x ^/ +3"^, and average 
if-band magnitudes are derived from sine-curve fitting. 

A close inspection of periods of LPVs lying away from 
the Mira PL relation indicated a problem of aliasing pe- 
riods of ~ 1/2 and ~ 1 year. For a total of 39 vari- 
ables the most significant period in Fourier periodogram 
was shorter, approximately 1/2 year, but they could be 
equally well fit by twice as long periods. These are listed 
by Rejkuba l|2()()3l Tab. 1). Other 7 variables have possible 
shorter periods and for 6 more LPVs improved amplitudes 
and mean magnitudes have been determined, but their pe- 
riods remained unchanged to within 3%. In the following 
analysis periods, amplitudes and mean magnitudes from 
the LPV catalogue (Rejkuba et al. I2()()3|l will be used ex- 
cept for those hsted in Tab. 1 of Rejkuba H2003|l . However, 
the results of our analysis are not influenced by this choice. 

3. LPVs in near-IR color-magnitude diagrams 

In Fig. n we show Kg- vs. Jg — Kg color-magnitude dia- 
grams (CMDs) for all the stars detected in at least 3 Kg- 
band epochs as well as in Jg and i/-band images that had 
ALLFRAME photometric errors smaller than 0.25 mag. 





Fig. 1. Top panels: Ks vs. Js — Ks color-magnitude diagrams for all the stars in Field 1 (left) and Field 2 (right) 
is shown with black dots. Bottom panels: Ks vs. H — Ks color-magnitude diagrams for Field 1 (left) and Field 2 
(right) stars. Variable stars with reliably determined periods (significance parameter from Fourier analysis < 0.7) are 
indicated with red circles and green boxes indicate those with less reliable periods (significance > 0.7). Dotted line at 
Ks — 21.24 shows RGB tip magnitude. 
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Field 1 CMD is displayed in the left and Field 2 in the 
right panel. Long period variable stars with reliably deter- 
mined periods (significance parameter from Fourier analy- 
sis < 0.7) are indicated with red circles and green boxes in- 
dicate those with less reliable periods (significance > 0.7). 
In the following analysis only those variables with bet- 
ter determined periods (significance < 0.7) will be used. 
These are 280 LPVs in Field 1 and 617 LPVs in Field 2. 

Dotted line at Kg — 21.24 shows RGB tip magnitude 
measured from the if-band luminosity function (Rejkuba 
I2003|l . There is a plume of foreground Galactic stars 
that are bluer and well separated from the red giants in 
NGG 5128. According to field simulations from Marigo 
et al. H2003|l they are mainly disk stars. In particular, old 
disk turn-off stars are expected to be distributed mainly at 
(Js — Ks) — 0.36, Galactic RGB and red clump stars have 
colors around {Jg — Kg) ~ 0.65 and low mass dwarfs with 
M < 0.6 Mq are found around ( Js - K^) ~ 0.85. The fact 
that they are distributed in vertical sequences is entirely 
due to a range of distances that they span. No Galactic 
foreground contamination is expected at (Js — Ks) > 1.0, 
where red giants in NGG 5128 are located. The equivalent 
Ks vs. H — Kg GMDs are shown in the two lower panels 
of Fig. n Here there is some overlap between the regions 
where foreground sources and red giants in NGG 5128 
are located. It should be noted from these GMDs that, 
due to the much higher density of stars in Field 2 and 
slightly worse average seeing during the observations and 
higher total luminosity sampled, the limiting magnitude 
is brighter and the photometric scatter is much larger. 

The i^s-band magnitude used to compose Jg — Kg and 
H — Kg colors in Fig. ^ is a single epoch measurement 
obtained closest in time to the respective Jg or _ff-band 
epoch. The time separation between the two is less than 
10 days. In contrast, the Kg magnitude plotted on the 
y-axis is a mean value from all the epochs and it well 
represents mean magnitudes of AGB stars. If a mean Kg 
magnitude is used for Jg ~ Kg and H ~ Kg color, a much 
larger scatter in color is observed in the part of the GMD 
where variable stars are found. 



100 



50 




-Fl 
F2 



0.5 



1 1.5 

Amplitude 



Fig. 2. Histogram showing the amplitude distribution of 
LPVs in Field 1 (sohd hne) and 2 (dotted line). 




Fig. 3. Amplitudes as a function of periods (top) and 
mean magnitudes (bottom) for Field 1 (left) and 2 (right). 



4. Amplitude distribution 

Amplitudes of the LPVs in NGG 5128 were determined 
as peak-to-peak magnitude differences of the best fitting 
sinusoid light-curve. Fig. El shows the amplitude distri- 
bution of all the LPVs with well determined periods. 
The majority of variables have if -band amplitudes of 
around 0.7 mag, with median value of 0.77 in both Fields. 
The JsT-band amplitudes for LPVs in LMG (e.g. Wood et 
al. I1983|l and in the Solar neighbourhood (e.g. Whitelock 
et al. I2000|l range between 0.1 and 1.2 mag. Only a few 
LPVs in the Sgr I field in the Bulge have if-band ampli- 
tudes larger than 1 mag with the largest amplitude of 2.5 
mag. 

The dependence of the amplitudes on the mean mag- 
nitudes and on the periods is shown in Fig. |3 As was also 
found for Galactic Mira-like LPVs (Whitelock et al. .2000,1 . 



there is a very weak correlation in the sense that longer 
period variables tend to have larger amplitudes. Similarly 
Glass et al. fl995) found that the large-amplitude stars 
have periods longer than 350 days in the Sgr I field in the 
Galactic bulge. 

There are only 32 (out of 897) variables with ampli- 
tudes smaller than AK < 0.4. Such small amplitude vari- 
ables usually tend to have fainter magnitudes. Our cat- 
alogue is much less complete for these variables. While 
some of them may have been detected as variable stars, 
no reliable periods could be measured for them. Semi- 
regular variables have smaller amplitudes making the de- 
tection of periodic variability even more difhcult with a 
small number of observations. Additionally, some of them 
might not even be detected as variable due to combination 
of small amplitudes and photometric uncertainty. They 
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Fig. 4. K-hand luminosity function for all variables with 
well established periods (Fourier analysis significance < 
0.7) in Field 1 (sohd hne) and Field 2 (dotted line). 
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Fig. 5. Comparison of the iiT-band luminosity functions of 
Field 1 in NGC 5128 halo (beaded blue histogram) with 
those of the Galactic bulge (solid black histogram) and 
of two fields in the bulge of M31 (dotted magenta his- 
tograms). All luminosity functions have been normalized 
in the range —6.5 < Mk < —5.5. The Galactic bulge data 
are from Zoccali et al. II2U(J3|I and the two M31 bulge lu- 
minosity functions are from Stephens et al. ^20031 fields 
F3 and F4). 

might comprise the large majority of the "non- variable" 
AGB stars detected in Field 1 (see also the next section). 

5. Luminosity function 

The luminosity functions for all the variables with well 
determined periods is shown in Fig. 0] with solid line for 
Field 1 and dotted line for Field 2 LPVs. The luminosity 
functions in the two fields are very similar. 



A detailed calculation of the expected number of LPVs 
in a ~ 10 Gyr old stellar population according to Renzini's 
(PggF) work yields 250-500 LPVs in Field 1 and between 
350 and 1500 LPVs in Field 2 (Rejkuba '5003. The large 
uncertainty of the expected number of LPVs in Field 2 
comes from highly uncertain total luminosity of this field 
due to higher number of blends and difficulty to correctly 
estimate the sky background. The numbers of expected 
LPVs are in very good agreement with the 437 and 709 
LPVs with measured periods in Fields 1 and 2, respec- 
tively. This would suggest that the majority of the LPVs 
belong to an old and metal-rich population similar to that 
of our Galactic bulge. 

Taking into account the foreground Galactic contami- 
nation and a probable number of blends, LPVs account for 
at least 26% and 70% of the extended giant branch popu- 
lation in Fields 1 and 2, respectively fReikuba 12003(1 . The 
completeness of the LPV catalogue in the two fields is 
rather similar as a function of period and amplitude and 
thus there seems to be an excess of non-variable bright 
AGB stars in Field 1. At least some 1150 stars brighter 
than the RGB tip cannot be accounted for by LPVs, fore- 
ground stars or blends in Field 1. In Field 2 this num- 
ber is an order of magnitude lower and can well be ac- 
counted by the incompleteness of the variable star cat- 
alogue. Basically all the constant luminosity AGB stars 
in Field 1 are fainter than Kg ^ 19.8. As mentioned 
above they could be semi-regular variables with ampli- 
tudes up to about AK ^ 0.3 mag, the limit at which 
the LPV catalogue completeness starts to drop rapidly 
from ~ 70%. Virtually no variables with AK ~ 0.1 or 
smaller are expected to be detected according to the com- 
pleteness simulations (Rejkuba et al. 12003(1 . In Field 2 
small amplitude variables could be hidden among blends. 
It should be noted that in the Galactic bulge and in old 
and metal-rich globular clusters, no constant luminosity 
stars are expected to be found above the RGB tip (Frogel 
& Whitelock Unna Class & Schultheis l^nH^ . New, de- 
tailed searches for variable stars in the LMC indicate that 
the large majority of the stars above the RGB tip, as well 
as many fainter than the RGB tip, are variable (Wood 
Ennrn Ita et aLEpi, Kiss & Bedding I2SSI). 

The positions of the brightest LPVs have been care- 
fully checked on the best seeing images. They show that 
the three brightest Field 1 variables are much brighter 
than expected, due to blending. One of them, with < 
K >~ 18.1 and Js — Kg ~ 1.75, is actually blended with 
a background galaxy and the very blue bright LPV at 
Js — Ks = 0.55 is contaminated by a bright neighbouring 
foreground star. The brightest 6 LPVs in Field 2 do not 
seem to suffer any blending. 

The comparison of the Field 1 K-hand luminosity func- 
tion with the luminosity functions of the Galactic bulge 
(ZoccaH et al. "lUU^ and of two bulge fields of the M31 
(fields F3 and F4 from Stephens et al. 2003 1 is shown in 
Fig. |5| All luminosity functions have been normalized in 
the range —6.5 < Mk < —5.5. Taking into account the 
blends, NGC 5128 Field 2 luminosity function is similar 
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Fig. 6. Histogram showing the period distribution of 
LPVs in Field 1 (soUd Une) and 2 (dotted Hne). 

to that of the Field 1 (see also Reikuba OIMIl . The tip 
of the AGB is around Mk = —7.5 in the Galactic bulge. 
However, Mira variables are known in Baade's Window 
and in Sag I fields that can reach Mk = —8 (Glass et 
al.Hnnni Shulthels & Glass IW)T)l . similar to the AGB tip 
observed in the M31. In NGC 5128 there are many more 
AGB stars and the tip of the AGB is more than half a 
magnitude brighter, clearly implying the presence of the 
intermediate age population in the halo. The tip of the 
AGB, measured from the average Ks-^and magnitude of 
the LPVs, is = 19.3. Rejkuba H2()()3|l has measured 
distance modulus of NGC 5128 of 27.92 ± 0.19 mag, im- 
plying that the brightest AGB variables in NGC 5128, 
which are not brightened by blending, have Mk = —8.65. 
Bolometric corrections were calculated from Montegriffo 
et al. (|T998) and the brightest AGB variables in Field 1 
have Mboi — —5.3, while there are few Miras in Field 2 
that reach reach M^oi — —5.7 mag. The brightest LPV 
in Field 2 has actually Mk = -9.78 and Mtoi = -6.43. 
Its red colors suggest that it is a long period variable. Its 
period is 696 days and it lies ~ 0.8 mag above the Mira 
PL relation fR.eikuba I2003|l . Its amplitude is 0.36 mag, 
rather low for its long period. Unless it is a blend of two 
perfectly aligned stars (hence not visible on 0'.'36 image), 
it is a good candidate for a hot-bottom burner (Blocker & 
Schonberner I1991|l . 

6. Period distribution of LPVs 

The distribution of periods for all the variables with good 
periods is shown in Fig. The distributions for Field 1 
and 2 are similar at the long period end, but there is a 
clearly larger percentage of stars with periods in the range 
of 250-300 days in Field 2. This is the same period range 
as the old and metal-rich LPVs in Galactic globular clus- 
ters. Since the completeness of the variable star catalogue 
in this period range is of the order of 80-90%, the larger 
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Fig. 7. Comparison of period distributions for Miras in 
globular clusters (Feast et al. l2()()2|l . in the Galactic bulge 
Sg rI field (Glass et al. ITTTO and in the LMC (Cioni et 
al. IM7T)l with NGC 5128 Miras in both fields studied in 
this work. Optically discovered Miras in the SgrI field in 
the Galactic bulge are shown as open histogram, while the 
10 IRAS sources identified in the same field are plotted 
with shaded histogram. Miras {AK > 0.4) in NGC 5128 
with Js — Kg colors redder than Jg — Kg > 1.6 are shown 
with shaded histogram. 

number of short period variables in Field 2 indicates dif- 
ferences in the stellar populations. 

The shortest periods are 160 and 155 days, respec- 
tively. The longest reliable periods are 1097 in Field 1 and 
1229 in Field 2. Due to the limited time interval spanned 
by the observations, amounting to 1197 days, most of the 
variables with well determined periods, hence observed 
through at least one to two periods, have periods shorter 
than 600-700 days. If there is a larger population of vari- 
ables with periods in excess of 700 days, only those 
with larger amplitudes would be detected. The steep drop 
in the period distribution present already around P=500 
days indicates that there is no significant population of 
such very long period variables. In the present distribu- 
tion, LPVs with periods longer than 500 days make 10% 
and 11% of the population in Fields 1 and 2, respectively. 
For periods longer than 550 days, the contribution to the 
distribution drops to 4% and 6%. The average period of 
897 LPVs with well determined periods is 395 days. For 
52 variables, periods were revised by Rejkuba (2003). If 
instead, for all the stars periods from the LPV catalogue 
(Rejkuba et al. 12003(1 are used, the average period is 388 
days. 

There are 58 variables in Field 1 and 25 in Field 2 
for which best fitting periods were shorter than 150 days. 
However, all of these period determinations are less se- 
cure. Some short period variables, belonging to Cepheid 
class are expected to be found in Field 1 where a conspic- 
uous recent star formation is observed (Mould et al. i2000l 
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Rejkuba et al. I2(JI)1|I . Their colors would be bluer than 
those of Miras and semi-regulars. Most of the blue vari- 
ables detected in our fields are unfortunately blended with 
foreground stars. 

The period distribution of LPVs in NGC 5128 can be 
compared with period distributions of LPVs in other stel- 
lar systems (see Fig. Ej). In old and metal-rich Galactic 
globular clusters, Miras span a rather narrow range of pe- 
riods 140 ^ P ^ 310 (Feast et aLEDH^I, while local Miras 
have periods ranging from ~ 100 to ~ 550 days (Whitelock 
et al. Unnni). According to Whitelock et al. I^TW^l there 
is no significant population of Miras with periods longer 
than ^ 700 days in the Galactic bulge. However, their re- 
sult is in disagreement with Harmon & Gilmore H1988|l 
and Blommaert et al. H1998|l . These last authors found 
some OH/IR variables with periods in excess of 2000 days 
in the Galactic bulge. Glass et al. H1995|l and Schultheis 
& Glass {20015 have analysed Sgr I field and Baade's win- 
dow in the Galactic bulge. The average period of 70 LPVs 
in Sgr I field is Pav = 333 day. In the Galactic Centre 
field the average period of LPVs is 427 days, while that of 
the OH/IR stars amounts to 524 days (Glass et al. 12001;) . 
Apparently, longer-period, hence younger, LPVs are con- 
centrated more towards the Galactic Centre, where a mix 
of young and old populations is observed due to a recent 
(^ 1 Gyr old) starburst. 

Glass et al. p995f) discussed the dependence of 
the type of LPVs discovered on the technique, in the 
sense that optical searches found principally shorter pe- 
riod Miras, while infrared techniques (IRAS and OH/IR 
sources) preferentially found longer period Miras. In our 
near-IR search, most of the short and longer period Miras 
should be discovered, with the exception of the most dust 
enshrouded OH/IR sources, if these were present. We find 
in our fields as well that the periods of the redder Miras 
are typically longer as shown by the shaded histogram in 
Fig. 13 Unfortunately due to rather low statistics it is not 
possible at this moment to make similar comparisons in 
other systems. 

In the LMC LPVs (including Mira and SR variables) 
span a range of periods from about 20 to ~ 1000 days (e.g. 
Wood '50011. Cioni et al. IMm Kiss & Bedding 2003). 

It has been now quite firmly established that the 
period of Mira variables represents a good indicator 
of the stellar population to which it belongs (Feast & 
Whitelock .1987,) . Longer period Miras are expected to 
have higher mass progenitors and therefore belong to a 
younger population (e.g. Iben & Rcnzini T983 Jura & 
Kleinmann 1992a,. 1992bi Kerschbaum & Hron ,.1992. Feast 
& Whitelock I^OOOIl . 

Unfortunately a precise empirical calibration of the 
initial mass or, equivalently, the age vs. period of Miras 
is lacking. It can be either inferred from the kinemati- 
cal properties of Miras in the solar neighbourhood (Feast 
119631 Feast & Whitelock 12000(1 . or from the observations 
of Miras in simple stellar populations like star clusters. 
According to Feast & Whitelock H1987(l , the kinematic age 
of a 400 day Mira is of the order of ~ 5 Gyr and an 850 



_ Fl 




H-Ks 

Fig. 8. Color-color diagrams for all the stars with 
ALLFRAME photometry errors smaller than 0.25. LPVs 
are indicated with red circles (secure periods) and green 
squares (less certain periods). Reddening vector for E(B- 
V)=0.5 is shown. Lines indicate intrinsic colors of Galactic 
giants (solid red line) and dwarfs (dashed blue line) from 
Bessell & Brett ^SB- 



day period Mira would be younger than about 4 Gyr. Feast 
1)1996(1 on the other hand explored the metallicity depen- 
dence of the mass-period relation for Miras and suggested 
that a 400 day Mira with a solar metallicity would have a 
mass of ^ 1.4 Mq. According to models (Iben & Renzini 
il983.. Vassiliadis & Woodil993j) a 400 day Mira would have 
initial mass of about ~ 1.0 Mq. However, it should be re- 
membered that the age of the AGB star with a 1.0 Mq 
progenitor will strongly depend on metallicity. The more 
metal-rich, the older the star, assuming a constant mass. 
For example, the turn-off age for a Z=0.004 star of 1.0 Mq 
is ^ 4.6 Gyr, while a same mass solar-metallicity turn- 
off star would be 7.7 Gyr old (Pietrinferni et al. 2003, in 
preparation). Turn-off age for a solar metallicity star of 
1.4 Mq is less than 3 Gyr. 
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Miras in the old and metal-rich globular clusters have 
ages of the order of 10-12 Gyr, with corresponding masses 
of the order of 0.6— O.TMq and periods around ~ 250 days. 
Another empirical point, calibrating the mass-period rela- 
tion is given by Nishida et al. |2000), who have determined 
periods of 3 Miras in the LMC star clusters with ages of 
1.6 to 2.0 Gyr. All three Miras have similar periods rang- 
ing from 491 to 528 days and they follow the Mira PL 
relations determined from shorter-period stars (Feast et 

NGC 5128 has a similar period distribution to Galactic 
bulge Miras, but there is a tail with some ~ 10% of the 
Miras having periods longer than 500 days. Unless this is 
purely due to age-metallicity-period degeneracy, we con- 
clude that this is an evidence for an intermediate-age AGB 
component in the NGC 5128 halo. 




7. Are there C-stars in NGC 5128? 

7.1. Color-color diagram 

Color-color diagrams for all the stars with ALLFRAME 
photometry errors smaller than 0.25 are shown in Fig. |S1 
Magnitudes have been de-reddened adopting a foreground 
reddening of E(B - V) = 0.11 and a CardeUi et al. [T^ 
reddening law. A reddening vector corresponding to E(B— 
V) = 0.5 is plotted. As for Fig. ^ single-epoch /-iT-band 
magnitudes closest in time to J and H-hand measure- 
ments are used to form colors. The solid (red) line is a 
Bessell & Brett (1988) fiducial to the Galactic giants and 
the dashed (blue) line is for dwarfs. Bessell & Brett fil- 
ters have been transformed to the photometric system of 
ISAAC (Chris Lidman, private communication). 

Carbon stars arc typically found among intermediate- 
age metal-poor populations. Single stars enter the AGB 
phase as oxygen-rich (0-rich) and they can be converted 
to carbon-rich (C-rich) after the products of nuclear burn- 
ing are brought up to the surface during thermal pulses. 
Due to the lower oxygen abundance of a metal-poor star, 
fewer C-atoms are needed to convert it into a C-rich star 
than is the case for a more metal-rich star. The ratio of 
the number of carbon- rich stars to oxygen-rich stars (C/M 
ratio) can therefore be used to derive metallicity gradi- 
ents (Cioni & Habing, I2003|l . Detecting carbon stars in 
NGC 5128 is important because their presence would be a 
definite proof of an intermediate-age and metal-poor com- 
ponent. Their origin might belong to a recently accreted 
metal-poor, LMC-like companion galaxy. 

Carbon stars in the LMC have 1.5 < J ~ Kg < 2.0. 
Stars redder than J—Kg > 2.0 are obscured and can either 
be 0-rich or C-rich (Cioni et al. IM7T]l . In the H - K 
vs. J — H color-color diagram, carbon stars are located 
redwards of H - K > O.A and J - H > 0.8. However, O- 
rich LPVs can also be found there. According to Bessell & 
Brett (1988) there should be no 0-rich LPVs redder than 
- X ^ 1.0 and J - ^ 1.4 (see their Fig. A3). There 
are 11 stars in Field 1 and 15 in Field 2 with H — K and 
J — H colors consistent with being carbon stars and all 




H-K 



Fig. 9. Bona fide C-star candidates in NGC 5128: LPVs 
with H - K ^ 1.0 and J - H Z 1.5 in NGC 5128 
are plotted as filled blue circles (Field 1) and filled red 
squares (Field 2). Obscured 0-rich (open circles) and C- 
rich (crosses) LPVs in the LMC (Whitelock et al. HUU^ 
are plotted for comparison. See text for the explanation 
of the model colors. 



but three of them have periods longer than 450 days and 
the average period of these stars is 526 days. Their mean 
amplitude is 1.2 mag. All of them have Jg — Kg > 2.5. Very 
red Jg — Kg colors and large amplitudes indicate that they 
probably have circumstellar shells. 

A comparison of these stars with obscured C-rich and 
0-rich AGB variables in the LMC (Whitelock et aL OUMjl 
is shown in Fig. El The LMC 0-rich AGB stars are plot- 
ted with open circles and the C-rich stars with crosses. 
NGC 5128 candidate C-rich LPVs are plotted as filled 
blue circles (Field 1) and filled red squares (Field 2). In 
order to facilitate the comparison, ISAAC data have been 
transformed into the SAAO photometric system. Cohen 
et al. (1981") have shown that the {H - K) and (J - H) 
colors of optically visible C-rich stars in the Milky Way 
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are correlated due to line blanketing by molecular absorp- 
tion bands (shown as dashed line labelled "MW C-rich"). 
The correlation for the optically visible C-rich stars in the 
LMC (dashed line labelled "LMC C-rich" ) was published 
by Costa & Frogel (|1995|l , while the correlations for the ob- 
scured C-rich (green solid line labelled "MW C-rich" ) and 
0-rich (magenta solid line labelled "MW 0-rich") Milky 

EMI 




Fig. 10. Period-color relations. For comparison a linear 
fit to the LMC period-color relations (Feast et al. I1989|l 
is shown (solid line in all three panels). In period-(J — K) 
diagram a linear fit to Sgr I Miras (Glass et al. I1995|l is 
shown with dashed line. 



Way stars comes from Guglielmo et al. H1993|l . All four cor- 
relation equations are taken from van Loon 1)1999(1 . who 
transformed them to {H — K) and (J — K) on the SAAO 
system. Stars lying along the correlations for obscured O- 
rich and C-rich stars could belong to either category, but 
the majority lies along the lines where optically visible 
C-stars are found. These are bona fide C-star candidates. 
However, more C-stars could be present among the long- 
period stars with somewhat bluer colors, but without spec- 
troscopy they cannot be distinguished from the obscured 
0-rich LPVs. 



7.2. Period-color distributions 

Color-log period diagrams for NGC 5128 LPVs arc shown 
in Fig. ^1 For comparison, a linear fit to the LMC color- 
log period relations (Feast et al. I1989|l are shown as solid 
lines. They have been shifted to the ISAAC photomet- 
ric system. It should be noted that these linear relations 
are valid only for P < 420 day. As for Figs. ^ and |H1 
single-epoch K-hand magnitudes closest in time to J and 
i?-band measurements were used to form colors. This pro- 
duces a large scatter around the mean Period-color rela- 
tion. The scatter is the largest in J — H vs. log(P) dia- 
gram where the phase difference between the observations 
is larger and thus some colors might even be non-physical. 

(J — H)o and (J — K)o colors are likely to be affected 
differently by metallicity changes, but also by gravity and 
atmospheric extension. Whitelock et al. H1991|l have com- 
pared (J — H)o and (J — K)o vs. logP diagrams for the 
LMC, the Galactic globular clusters, the solar neighbour- 
hood and the Galactic bulge. Their conclusion is that at 
a given period, colors of Miras in these different environ- 
ments are very similar. However, Miras in the Sgr I field of 
the Galactic bulge (Glass et al. ll995|l have redder {H—K)o 
and slightly bluer (J — H)o colors than the LMC Miras at 
a given period, while (J — K)q colors do not show so large 
offset. Feast ()199f)|l reports the mean {J—K)o color-period 
relation for Sgr I Miras studied by Glass et al. (|1995)I 
which has quite a different slope with respect to the LMC 
relation. He ascribes the difference to metallicity differ- 
ences. The SgrI relation is shown in the (J — K)o color- 
log P diagram as a long dashed line. It provides much bet- 
ter fit to NGC 5128 LPVs and suggests a similar metallic- 
ity as in the Galactic bulge. The large scatter around the 
color-period relation due to random phase J^-band obser- 
vations and possible abundance spread within NGC 5128 
halo (e.g. Harris et al. 119991 Rejkuba et al. 20^ prevents 
us from drawing more quantitative comparison. 
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Some very red stars are present among the longest pe- 
riod variables in NGC 5128. This has also been seen in 
the LMC, the solar neighbourhood and the Bulge (Glass 
et al. I1995L Kiss & Bedding i2003)) . and is ascribed to the 
presence of cool circumstellar shells arising from mass loss. 
Such shells contribute principally at longer wavelengths 
creating an infrared excess. Actually there are 113 sources 
in Field 1 and 236 in Field 2 with Kg-hand magnitudes 
brighter than the RGB tip and no Js or iJ-band coun- 
terparts. Many of these could be dust enshrouded LPVs 
similar to those found in the central parts of the Milky 
Way by van Loon et al. H2Ul)3|l . Unfortunately, periods 
could only be measured for a handful of them. 

8. Conclusions 

We have analysed near-IR properties of 897 LPVs with 
well determined periods and JgHKs photometry in two 
halo fields in NGC 5128. Mostly they are found to be 
brighter than the tip of the RGB {Kg < 21.24; Rejkuba 
I2003|l with magnitudes ranging from about ii' = 19 to 
K = 21.5. They have periods between 155 and 1000 
days and K-band amplitudes between 0.1 and 2 mag, 
characteristic of semi-regular and Mira variables. They 
obey the same period-luminosity relation as found in the 
Magellanic Clouds, Solar neighbourhood and the Galactic 
bulge (Rejkuba '2003). Here we have compared the col- 
ors, periods and amplitudes of these variables with those 
found in old stellar populations like Galactic globular clus- 
ters and Galactic bulge as well as with intermediate-age 
Magellanic Cloud long period variables. 

There is an excess of short period, globular cluster-like, 
Mira variables in Field 2 showing the differences between 
the stellar populations of the two fields. The brightness 
of the AGB tip, which occurs around Mk = —8.65 in 
NGC 5128, corresponding to Ahoi — —5.3, compared to 
AIk = —8 in the Galactic bulge (Zoccali et al. l2003;l and in 
the M31 bulge, (Stephens et al. l2003|l is an evidence for a 
~ 4 Gyr old component (Mould & Aaronson ll982(l . Some 
10% of the long period P > 500 day Miras are expected to 
have ages younger than 7 Gyr and could be as young as 3— 
4 Gyr old. Approximately two dozen C-star candidates are 
identified, but more could be present among the reddened 
long period variables. C-stars are expected to be found 
only among metal-poor AGB stars. 
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